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Abstract 



We report on a study of the invariant mass spectrum of the hadronic system 
in the decay r~ — > 7r~7r*^z/,-. This study was performed with data obtained 
with the CLEO II detector operating at the CESR e^e~ colhder. We present 
vQ I fits to phenomenological models in which resonance parameters associated 

■^ ■ with the p(770) and p(1450) mesons are determined. The tt^tt^ spectral 

^^ ■ function inferred from the invariant mass spectrum is compared with data on 

e^e~ — > 7r+7r~ as a test of the Conserved Vector Current theorem. We also 
^^ I discuss the implications of our data with regard to estimates of the hadronic 

■"^ ' contribution to the muon anomalous magnetic moment. 
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I. INTRODUCTION 

The r is the only lepton heavy enough to decay to final states containing hadrons. Since 
leptons do not participate in the strong interaction, r lepton decay is well-suited for isolating the 
properties of hadronic systems produced via the hadronic weak current [|^J2|. Furthermore, angular 
momentum conservation plus the transformation properties under parity and G-parity of the vector 
and axial vector parts of the weak current give rise to selection rules that constrain the types of 
hadronic states that may form. Thus, r lepton decay provides an especially clean environment for 
studying these states. In this article, we present a study of the 7r~7r'' system produced in the 
decay r~ — > 7r~7r°i/^ based on data collected with the CLEG II detector. 

In semi-hadronic r decay, hadronic states consisting of two pseudoscalar mesons may only have 
spin-parity quantum numbers J^ = 0"^ or 1~. In addition, the Conserved Vector Gurrent theorem 
(GVG) forbids production of 0^ non-strange states in r decay. Thus, within the picture of resonance 
dominance in the accessible range of squared momentum transfer g^, the decay r^ -^ tt'tt^Ut is 
expected to be dominated by production of the lowest lying vector meson, the p{770). Radial 
excitations, such as the /o(1450) and the />(1700), may also contribute. Although these are well- 
known mesons, their properties have not been measured precisely, and there exists a wide variety of 
models that purport to characterize their line shapes. New data can help improve the understanding 
of these states. 

Finally, GVG relates properties of the vr^vr^ system produced in r decay to those of the vr+vr" 
system produced in the reaction e^e~ -^ vt+tt^ in the limit of exact isospin symmetry. The degree 
to which these relations hold has important consequences. For example, data on the e~^e~ process 
is used to determine the dominant contribution to the large but uncalculable hadronic vacuum- 
polarization radiative corrections to the muon anomalous magnetic moment a^ = {g^ — 2)/2. With 
GVG, T data can be used to augment the e^e~ data, leading to a more precise Standard Model 
prediction for the value of a^ [Q . 

Here, we attempt to address some of these issues, using a high-statistics, high-purity sample of 
reconstructed t~ — > 'K~T:^Vr decays. The measurements presented here supercede earlier prelimi- 
nary results from GLEG II on this subject 0. Work in this area has also been published by the 
ALEPH Gollaboration H. In Sec. |ll|, we review models of the hadronic current in the decays of the 



r to vector mesons, and specify the models we employ to extract resonance parameters. In Sec. Ill 



we discuss our data sample and the event selection criteria. To mitigate experimental biases, we 



apply several corrections to the data, described in Sec. IV. The results of fits to the corrected 



q^ spectrum are reported in Sec. ^, and systematic errors are discussed in Sec. VI. We compare 



our data with those obtained by ALEPH and the low-energy e^e experiments in Sec. VII. In 



Sec. VIII we discuss the applicability of our data for predictions of the muon anomalous magnetic 



moment. Finally, we summarize our results in Sec. IX 



II. PHENOMENOLOGY AND MODELS 
A. Model-Independent Phenomenology 

The decay rate for r~ -^ ■k~'k^i't- can be written as ||l| 



dr{r- ^ 7T-7T^l^r) _ G^lK^r^^V ,,^2 2^2 



dg2 327r2 M^ 
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{M^-qy{M^ + 2q')v-- {q'), (1) 



where q^ is the invariant mass squared of the 'k~-k^ system, and v'"'^ (q^) is the vector spectral 
function characterizing the (a priori unknown) hadronic physics involved in the formation of the 
(J = 1~) 7r~7r'^ system. Gp is the Fermi constant, Vud the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element, and Mj- the r lepton mass. S^^ denotes electroweak radiative corrections not 
already absorbed into the definition of Gp, some components of which have been determined 
theoretically 0-^. 

The corresponding vr+vr" spectral function v'^'^{q^) can be inferred from the cross section for 
e+e" ^vr+vr" ^^■. 

/4 2 2 \ 
a(e+e"^7r+^-)=P^^^U""(s), (2) 

where s = cf' Is the squared e^e~ center-of-mass energy. Up to isospin- violating effects, CVC 
allows one to relate the spectral function obtained from r decay to the isovector part of the e'^e~ 
spectral function: 

vr=Aq^) = v--\q') . (3) 

The e^e~ spectral function can also be expressed in terms of the pion electromagnetic form 
factor F^(g^): 



1 



2p7r 
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-^^{^') = TKZ FM') ^ h (4) 
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where the last factor represents the P-wave phase space factor, with Pt^ being the momentum of 
one of the pions in the tttt rest frame. The r decay spectral function can be similarly expressed in 
terms of the weak pion form factor. 

B. Models of the Hadronic Current 

The hadronic physics is contained within v^'^ic^^, or equivalently Fj^^q^). From the electric 
charge of the vr", it is known that -F7r(0) = 1. Beyond that, its form at low energies is not presently 
calculable in QCD, and models must be used. With resonance dominance, it is expected that F^ 
is dominated by the line shape of the /o(770) meson, with contributions from its radial excitations, 
the /9(1450) and /9(1700) mesons (denoted as p' and p" , respectively). 

Various Breit-Wigner forms have been proposed ||l|Jll|-14| to parameterize i^V- We consider 



here two models: those of Kiihn and Santamaria |l^ and Gounaris and Sakurai [11|, denoted as 
the K&S and G&S models, respectively. 



1. The Model of Kiihn and Santamaria 
In addition to its simplicity, the K&S model is useful since it is implemented in the TAUOLA r 



decay package [15| used in the CLEO II Monte Carlo simulation. The form is given by 



Py^'^il') = i^f^^^ + ... (^^P + PBWp, +^BWp. + •••), (5) 



where 



BWr, 



M2 



(M2 



/?r,(g2) 



(6) 



represents the Breit-Wigner function associated with the yo(770) resonance hne shape, with Afp and 
Tp{q^) denoting the p meson mass and mass-dependent total decay width. The assumed form for 
the latter is described below. The parameters (3 and 7 specify the relative couplings to p' and p" , 
and the ellipsis indicates the possibility of additional contributions. The Breit-Wigner functions 
are individually normalized so that the condition -F7r(0) = 1 is satisfied with the inclusion of the 
1/(1 -|- /? + 7) factor. For application to e^e~ -^ vr+vr^ data, one must consider isoscalar as well 
as isovector contributions. For this, the form for FtI ~ ' ' is obtained by modifying BWp so as to 
characterize p-co interference, which is not relevant for r decay. 

An alternate form for v^'^{q'^) can be obtained from consideration of the amplitudes for weak 
production and strong decay of p mesons. For the case where only the p{770) contributes, the 
spectral function can be expressed (following Tsai M) as: 



z;""(g2 



2vr/2 



'rp{q')/n 



(M2-g2)2 + ^2r2(g2) 



where 



^^(^^)- US 




(7) 



(8) 



gives the energy dependence of the p width. The constants fp (with units of mass squared) and 
Qp (dimensionless) can be identified as the weak and strong p meson decay constants, respectively. 
The -F7r(0) = 1 condition is satisfied for fpQp = \/2Mp, in which case the K&S form is recovered. 

The energy dependence of the p width may be more complicated than the P-wave behavior 
indicated in Eq. g. Various authors [12,16| suggest the need for an additional Blatt-Weisskopf 
centrifugal barrier factor |17| which takes the form 



2^2 



'R 



1 + R'p'o 

1 + R^pI' 



(9) 



where po = Ptt{(1^ = Mp), and R denotes the range parameter with a value assumed to be of 
0(1 fermi/?ic). This factor multiplies the right hand side of Eq. pi and thus modifies the T{q'^) 
factors appearing in both the numerator and denominator of the Breit-Wigner form for v^'^{q'^) 
given by Eq. |^. 



2. The Model of Gounaris and Sakurai 



The G&S model |11| has been used by a number of authors [|, 11, 13, 18 1 to parameterize the 
e^e~ — > vr+TT^ cross section. In this model, the form for Fj^ is derived from an assumed effective 
range formula for the P-wave vr-vr scattering phase shift, assuming p(770) meson dominance. This 
yields 



FAq') 



MJ + dMpTp 



{M^-q2)+f^q2 



q^Tpiq^ 



(10) 



where Fp denotes rp(g^ = M^), and 



2\ i,/;\,f2> 



f{g')=pi{q') h{q')-h{M'p) -p^^(q^-M'p) 



^2 /'„2 )\,f2> 



dh 



P' rlr,2 



dq^ 



(11) 



%^)^^^ln^^^±#^, (12) 



pI tt^ 2M,, 

and d is chosen so as to satisfy the -F7r(0) = 1 condition, 

.= Mfin^P + ^^o+J^_J^. (13) 

7rp2 2AU 27rpo vrpj] ^ ^ 

Following Refs. ||6|, p!3| , p!^ ] , we employ an extension of this model to include possible p' and p" 
contributions, as in Eq. |5|. 

The G&S form for F-,^ is similar to the K&S form in that (1) both are normalized so that 
F^(0) = 1, and (2) their shapes are similar in the vicinity of the p peak, since /{q"^) in Eq. ^ goes 
as Mp — q^ near q'^ = M^ [^. However, the additional term in the numerator of Eq. |l^ results in 
a larger value for Fj^ at q"^ = M^ relative to that in the K&S model, given the same values for Mp 
and Tp. For Mp = 0.775 GeV, the value of d is 0.48, such that F^{M^) is larger by 9% than the 
corresponding value from the K&S model. 



III. DATA SAMPLE AND EVENT SELECTION 
A. Detector and Data Set 

The analysis described here is based on 3.5 fb^^ of e^e^ collision data collected at center-of- 
mass energies 2£'beam of ~ 10.6 GeV, corresponding to 3.2 x 10^ interactions of the type e"'"e~ ^• 
T"^r~(7). These data were recorded at the Cornell Electron Storage Ring (CESR) with the CLEO II 
detector [^9| between 1990 and 1994. Charged particle tracking in CLEO II consists of a cylindrical 
six-layer straw tube array surrounding a beam pipe of radius 3.2 cm that encloses the e~^e~ collision 
region, followed by two co-axial cylindrical drift chambers of 10 and 51 sense wire layers respectively. 
Scintillation counters used for triggering and time-of-flight measurements surround the tracking 
chambers. For electromagnetic calorimetry, 7800 CsI(Tl) crystals are arrayed in projective and axial 
geometries in barrel and end cap sections, respectively. The barrel crystals present 16 radiation 
lengths to photons originating from the interaction point. 

Identification of t~ — > 7r~7r''i'T- decays relies heavily on the segmentation and energy resolution 
of the calorimeter for reconstruction of the vr". The central portion of the barrel calorimeter 
(I cos^l < 0.71, where 9 is the polar angle relative to the beam axis) achieves energy and angular 
resolutions of aE/E{%) = 0.35/£^°-^5 + 1.9 - O.IE and cj^ (mrad) = 2.8/ VE + 2.5, with E in 
GeV, for electromagnetic showers. The angular resolution ensures that the two clusters of energy 
deposited by the photons from a vr'^ decay are resolved over most of the range of vr'^ energies typical 
of the r decay mode studied here. 

The detector elements described above are immersed in a 1.5 Tesla magnetic field provided by 
a superconducting solenoid surrounding the calorimeter. Muon identification is accomplished with 
plastic streamer tubes, operated in proportional mode, embedded in the flux return steel at depths 
corresponding to 3, 5 and 7 interaction lengths of total material penetration at normal incidence. 



B. Monte Carlo samples 

We have generated large samples of Monte Carlo (MC) events for use in this analysis. The 



physics of the r-pair production and decay is modelled by the KORALB/TAUOLA event generator |15|, 
while the detector response is handled with a GEANT-based [^^ simulation of the CLEO II detector. 
The primary MC sample, denoted as the generic r MC sample, consists of 11.9 million r-pair 
events with all decay modes present. We generated an additional sample enriched in r~ -^ tt~tt^i't 
decays, bringing the total number of MC signal decays to 10.9 million, corresponding to roughly 
seven times the integrated luminosity of the data. The generic r Monte Carlo sample is used 
to estimate backgrounds from non-signal r decays, as well as for comparisons of kinematic and 
detector-related distributions with those from the data. We employ the full MC sample for the bin 



migration and acceptance corrections described in Sec. IV. 

The 7r~7r'^ spectral function implemented in the Monte Carlo is the K&S model, with parameters 
{Mp, Tp, P, Mpi, Vp,) = (0.773, 0.145, -0.145, 1.370, 0.510) in GeV, except for (3 which is 
dimensionless. Here Tp denotes the pole p meson width, Tp[q^ = M^). These parameters are based 



on one of the fits by Kiihn and Santamaria [13| to the e'^e -^ vr+vr data. This fit did not allow 
for a possible p(1700) contribution. 

C. Event selection 

Tau leptons are produced in pairs in e^e^ collisions. At CESR beam energies, the decay 
products of the r+ and t~ are well separated in the CLEO detector. The decay of the r~ lepton 
into TT'ir^UT is referred to as the signal decay, while that of the recoiling r"^ is referred to as the 
tag decay, and similarly for the charge conjugate case. Due to limited charged tt/K separation 
capabilities, we do not attempt to distinguish tt" from K~ in this analysis. As a result, our selected 
event sample contains background from the Cabibbo-suppressed channel r~ -^ K~TT^Vr- This and 
misidentified decays from other channels are subtracted statistically using the generic r Monte 
Carlo sample described above. 

To reject background from non-rr events, we require the tag decay products to be identified 
with one of three decay channels: e'^Vei'T ("e tag"), fi'^UpV-r ("/i tag"), and 7r"'"7r%^ ("p tag"). For 
the "/9 vs. p" topology, each event is considered twice, corresponding to the two ways of labelling 
the decays as tag and signal decays. Thus, in such events both decays are used in our analysis 
if the requirements given below are met for both combinations of tag and signal labels. We have 
previously used these event topologies to measure the branching fraction for the signal decay mode, 
described in Ref. |^l|. The event selection used here is similar and is described below. 

We require an event to contain exactly two reconstructed charged tracks, separated in angle by 
at least 90°. Both tracks must lie in the central region of the detector: the tag track must lie within 
I cos^l < 0.8, while the signal track must have | cos^| < 0.71, so as to avoid excess interactions in 
the main drift chamber end plate. Both tracks must be consistent with originating from the e^e~ 
interaction region, and have momentum between O.OS-Ebeam and 0.90£^beam- The momenta of all 
charged tracks are corrected for dE/dx energy loss in the beam pipe and tracking system. 

Clusters of energy deposition in the calorimeter are considered as candidates for photons from 
vr'^ decay if they are observed in the central part of the detector (| cos9\ < 0.71), are not matched to 
a charged track, and have energy greater than 50 MeV. Pairs of photons with invariant mass M^^ 
within 7.5 a-y-y of the vr*^ mass are considered as vr*^ candidates. The 77 invariant mass resolution a-y^ 

8 



varies from 4 to 7 MeV/c^, depending on vr'^ energy and decay angle. The vr^ energy is required 
to be greater than O.OSE'beam- Each vr'^ candidate is associated with the charged track nearest 
in angle to form a ir~7r^ candidate. If more than one vr'^ candidate can be assigned to a given 
track, only one combination is chosen, namely that for which the largest unused barrel photon-like 
cluster in the 7r~7r'' hemisphere has the least energy. A cluster is defined to be photon- like if it has 
a transverse energy profile consistent with expectations for a photon, and if it lies at least 30 cm 
away from the nearest track projection. 

As mentioned earlier, backgrounds from multihadronic {e~^e~ — > qq) events are rejected by 
identifying the tag system as being consistent with r"*" decay to neutrino(s) plus e"*", //"•" or tt+tt^. 
The tag track is identified as an electron if its calorimeter energy to track momentum ratio satisfies 
0.85 < E/p < 1.1 and if its specific ionization in the main drift chamber is no more than two 
standard deviations (a) below the value expected for electrons. It is classified as a muon if the 
track has penetrated at least the innermost layer of muon chambers at 3 interaction lengths. If the 
tag track is not identified as an e or a //, but is accompanied by a second tt^ of energy > 350 Me^, 
then the track- vr*^ combination is classified as a p tag. The invariant mass of this combination must 
be between 0.55 and 1.20 GeV. 

To ensure that these classifications are consistent with expectations from r decay, events are 
vetoed if any unused photon-like cluster with | cos9\ < 0.95 has energy greater than 100 MeV, or if 
any unmatched non-photon-like cluster has energy above 500 MeV. Finally, the missing momentum 
as determined using the vr^vr'^ and tagging systems must point into a high-acceptance region of 
the detector (| cos^missl < 0.85), and must have a component transverse to the beam of at least 
0.08£^beam- These requirements also limit the misidentification of r decays containing multiple vr'^'s 
as signal decays. 

D. Final event sample 

With this selection, 103522 events remain. The distribution in normalized di-photon invariant 
mass 5^^ = (M^^ — M^o)/cr^^ for these events is shown in Fig. ||, with the corresponding Monte 
Carlo distribution overlaid. Of these, 94948 lie in the vr*^ signal region, defined to be the interval 
—3.0 < S^^ < 2.0. The asymmetry of the distribution and the signal region definition arises 
because of the asymmetric energy response of the calorimeter. The low-side tail of the photon 
energy response curve is due primarily to rear and transverse leakage of high energy showers out 
of the Csl crystals whose energy depositions are summed in determining the energy of a given 
photon. We also make use of 2281 events lying in the side-band regions —7.5 < S^^ < —5.0 
and 3.0 < S-y^ < 5.5 to model backgrounds associated with spurious vr" candidates. After these 
selections, we redetermine the photon energies and angles making use of the vr*^ mass constraint, 
so as to improve the iT~7r^ invariant mass resolution. 

The M^-^o spectrum is shown, after side-band subtraction, in Fig. y. The agreement between 
data and MC spectra is more than an indication of the validity of the application of CVC. It 
also suggests that the event kinematics in the MC samples are sufficiently similar to those in the 
data that the selection criteria described above are not likely to have introduced significant biases. 
Additional support for this is the comparison between data and generic Monte Carlo samples of 
the 7r~ momentum and vr" energy distributions, shown in Figs. H and HI Several events in Fig. lie 
above the r lepton mass. The small number of these events indicates that possible backgrounds at 
high mass, such as low-multiplicity qq events, are not significant. 
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FIG. 1. Normalized M^^ distribution for vr*^ candidate photon pairs in the data (points) and 
the generic r Monte Carlo sample (line histogram), after all other cuts. 
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FIG. 2. Raw M^-^o spectrum after vr" side-band subtraction for candidate decays from the 
data (points) and the generic r Monte Carlo sample (line histogram). 
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FIG. 3. Distribution in the momentum of the tt~ , divided by the beam energy, for candidate 
decays, from the data (points) and the generic r Monte Carlo (Une histogram) samples, after ir^ 
side-band subtraction. The bottom plot gives the deviations of the data spectrum from the Monte 
Carlo spectrum, normalized by the Monte Carlo spectrum. 
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After 7r° side-band subtraction, backgrounds from non-signal r decays are estimated to be 
6.63 lb 0.20% from the generic r Monte Carlo sample. The dominant channels are t~ — > t:~t:^'k^i't 
(4.01 ± 0.08%), T- -^ K'Ti^Vr (1.86 ± 0.16%) and r" -^ h-Rlir^i^r (0.59 ± 0.08%), where h 
denotes tt or i^ and the errors include branching fraction uncertainties as well as statistical errors. 

IV. CORRECTIONS TO THE DATA 

One goal of this analysis is to analyze the 7r~7r'' mass spectrum in the context of several models. 
However, we are not able to explore all possible models. In addition, it is desirable to compare our 
spectrum to data from other experiments in a model independent way, as well as to present it in a 
form that facilitates comparison with future data or models. These considerations motivate us to 
construct a histogram of the mass spectrum that has been corrected for (primarily) experimental 
effects. We then carry out simple x^ fits to the corrected spectrum using the models of the tt~'7t^ 
line shape described earlier. 

Three experimental effects give rise to distortions in the tt~tt^ mass spectrum: (1) backgrounds; 
(2) smearing due to resolution and radiative effects; and (3) mass-dependence of the experimental 
acceptance. In this section, we describe the corrections, that we applied in the order listed to 
remove these distortions. These corrections rely on the Monte Carlo simulation of the physics and 
detector response. 

A. Binning of the M^-^o Spectrum 

Before discussing the corrections mentioned above, we note that we have elected to bin the 
M^-j^o spectrum in intervals of 25 MeV below 1 GeV, and 50 MeV above 1 GeV. This binning is 
chosen so as to be sensitive to rapidly varying regions of the spectrum while limiting the size of 
the bin migration correction and consequently the magnitude of correlations among nearby bins in 
the corrected spectrum. This is important for the stability and accuracy of the x^ fit procedure, 



which is known to be biased when data points are strongly correlated [26|. The increase in mass 
resolution from approximately 6 MeV at low masses to 17 MeV at high masses motivates the large 
bin width above 1 GeV. The large bin width is also beneficial in the very high mass bins where 
low statistics could lead to non-Gaussian fiuctuations. 



B. Corrections for Backgrounds 

As noted earlier, the backgrounds entering the r~ — > -k'ti^Vt sample are small. Side bands in 
the M^^ distribution are used to model the fake-vr*^ contribution. The remaining 6.6% are modeled 
with the generic r-pair Monte Carlo sample. These subtractions are performed bin by bin in the 
mass spectrum. The Monte Carlo spectra for the signal and primary background modes are plotted 
in Fig. ^, after vr" side-band subtraction. 

C. Correction for Bin Migration 

Detector resolution causes the ■k~t[^ mass spectrum to become broader. The presence of radi- 
ation in the decay r~ -^ n^ -k^ {'^)vt- is also important. The radiative photons tend to be low in 
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FIG. 5. Contributions to the reconstructed M^~^o spectrum from Monte Carlo simulation of 
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FIG. 6. Effect of bin migration correction: M^-^o spectrum prior to (solid histogram with 
points), and after (dashed histogram) unfolding of resolution and radiative distortions. 



energy, and are difficult to distinguish from photons from a possible second n^ or from fragments of 
the hadronic shower from charged pions that interact in the calorimeter. Consequently, we can not 
identify them reliably, either for inclusion in the invariant mass calculation or as a basis for vetoing 
events. The net effect of ignoring decay radiation is to broaden and shift the mass spectrum. 

We correct for these effects by performing an approximate unfolding procedure, based on a bin 
migration matrix determined from the r~ — > 7r~7r'^^'T- MC sample. This procedure is outlined in 
Appendix^. Since the experimental resolution on M^^o is dominated by that on the vr'^ energy, the 
agreement between data and MC shown in Fig. [l| gives us confidence in this aspect of the correction 
procedure. For the radiative effect, we rely on the PHOTOS-based simulation |22| employed by 
TAUOLA. The unfolded spectrum is shown as the dashed histogram in Fig. p, with the uncorrected 
spectrum overlaid. 
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FIG. 7. Acceptance as a function of generated 7r~7r'^(7) mass, as determined from the full 
r^ -^ tt^tt^i't- Monte Carlo sample. 

D. Correction for Acceptance 

Finally we correct for mass dependence of the acceptance, plotted in Fig. |^. Again, this is 
determined from the MC simulation. The main effects causing this dependence are associated with 
the kinematics of the decay and the cuts imposed in the event selection, both of which are well 
understood. 

E. The Corrected Mass Spectrum 



The corrected mass spectrum is given in tabular form in Appendix O, along with elements of the 
covariance matrix characterizing the statistical errors and the correlations among entries introduced 



by the bin migration correction procedure. The spectrum is also made available electronically |23| 



V. RESULTS OF FITS FOR RESONANCE PARAMETERS 
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A. Fitting Procedure 

We perform ^ fits to the fully-corrected 7r~7r'^ mass spectrum to extract resonance parameters 
and couplings. The x^ minimization and parameter error determination is carried out using the 
MINUIT program [gj]. Because of poor statistics and/or uncertainties associated with the back- 
ground estimation and acceptance correction, only data in the range 0.5 to 1.5 GeV are included 
in the fits. Also, as a result of the unfolding procedure, off-diagonal terms of the covariance matrix 
are non-zero, and the corresponding terms must be included in the calculation of the x^ [^- The 
off-diagonal elements of the covariance matrix are not reflected in the error bars shown in the 
figures in this section. 

Since the functional forms used to fit the data are nonlinear functions of many parameters, 
several iterations of minimization are performed before convergence is reached. We also integrate 
the fit function within each bin when computing the x^- We have tested this procedure using 
high-statistics generator-level Monte Carlo samples to ensure the reproducibility and accuracy in 
the determination of fit parameters and their errors. 

B. Fit to the K&S Model 

In this section, we report in detail the results from the simplest fit, done using the K&S model 
with no p" contribution. Although this model is normalized according to the -F7r(0) = 1 constraint, 
we introduce an additional parameter multiplying the K&S function (Eq. S), which is allowed to 
float in the fit. We have elected to do this for several reasons. First, we do not expect this or 
any other model of the line shape of a broad resonance to hold arbitrarily far from its peak. If 
the model does not hold at very low values of M^-^o then enforcing the -F^(O) = 1 condition 
can bias the resonance parameters. Second, the focus of this analysis has been on the shape 
of the mass spectrum: for example, tight cuts have been applied to maintain high purity. The 



normalization of this spectrum (see Sec. VII B ) depends on external measurements which have 



experimental uncertainties, as well as on theoretical factors which also have uncertainties. Given 
these considerations, the value of strictly enforcing the normalization condition is questionable. 

The fully corrected M^^a spectrum with the fit function superimposed is displayed in Fig. @. 
The x^ for this fit is 27.0 for 24 degrees of freedom. We obtain 

Mp = 774.9 ± 0.5 ± 0.9 MeV, 
Fp = 149.0 ± 1.1 ± 0.7 MeV, 
/3 = -0.108 ±0.007 ±0.005, 
Mp, = 1364 ± 7 ± 8 MeV, 
Fp/ = 400 ± 26 ± 23 MeV, 

where the first error is statistical and the second is due to systematic uncertainties, described in 



Sec. Vl. When interpreted in terms of the pion form factor, the normalization gives |i^7r(0)|^ = 
1.16 ± 0.02, where the error is statistical only. Performing the same fit, but with the -FyrCO) = 1 
normalization condition imposed, yields a x^ of 62.1 for 25 degrees of freedom and significantly 
different values for the other fit parameters (i.e., Mp = 772.3 MeV, Tp = 144.6 MeV). 
The fit parameters are correlated, with the correlation matrix: 
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FIG. 8. Fully corrected M^-^o distribution in r~ -^ -k^ti^Ut events (points). The solid curve 
overlaid represents the results of the fit to the K&S model. The dashed curve is obtained using 
the p(770) parameters obtained from this fit, but with the p' contribution turned off (i.e., (3 set to 
zero) . 



TABLE I. Results from fits to the corrected M^-^o spectrum 


over the range 0.5-1.5 GeV, 


for several models. The errors 


shown are statistical only. K&S refers to the model of Kiihn 


and Santamaria 


13 while G&S refers to the model of Gounaris and Sakurai |Tl]. See text for 


descriptions of these models and their variants. 






Fit Parameter 






Model 






K&S 


K&S w/p" 


K&S w/barrier 


G&S G&S w/p" 


Mp (MeV) 


774.9 ± 0.5 


774.6 ± 0.6 


769.7 ±0.7 775.3 ±0.5 775.1 ±0.6 


Tp (MeV) 


149.0 ±1.1 


149.0 ±1.2 


145.8 ±1.3 150.5 ±1.1 150.4 ±1.2 


/5 


0.108 ±0.007 


-0.167 ±0.008 


-0.160 ±0.008 -0.084 ±0.006 -0.121 ±0.009 


Mp, (MeV) 


1364 ± 7 


1408 ± 12 


1321 ±9 


1365 ± 7 1406 ± 13 


Vp, (MeV) 


400 ± 26 


502 ± 32 


397 ±17 


356 ± 26 455 ± 34 


Rp (GeV-i) 


— 


— 


1.9 ±0.3 


— — 


Rp, (GeV-i) 


— 


— 


5.0 ±2.2 


— — 


7 


= 


0.050 ±0.010 


= 


= 0.032 ± 0.009 


\F.m'' 


1.16 ±0.02 


1.14 ±0.02 


0.39 ±0.01 1.04 ±0.02 1.03 ± 0.02 


xVdof 


27.0/24 


23.2/23 


22.9/22 


26.8/24 22.9/23 






/ 1.00 




\ 








0.64 


1.00 










0.75 


0.51 1.00 




(14) 






0.89 


0.68 0.43 


1.00 






0.21 


0.22 0.21 


0.16 1.00 








I -0.52 


-0.38 -0.08 - 


-0.74 0.33 1.00 j 





where the parameters are normalization, Mp, Tp, /?, Mp/, and Fp/, respectively. 



C. Fits to Other Models 



The results from fits of the corrected M^-^o spectrum to various models are given in Table |. 
Several of these fits are illustrated in Fig. p. For fits including a p" contribution, we fix its 
parameters to world average values [^] {Mpn = 1.700 GeV, T pn = 0.235 GeV), but allow the 
relative coupling constant 7 to float. 

The G&S fits behave similarly to the K&S fits, as suggested by the x^ values in Table ||, and 
by the nearly overlapping solid and dashed curves in Fig. ^. In the figure, the deviation between 
the K&S and G&S curves is only visible at very low and very high values of M^-^o. The deviation 
at low values is reflected by the difference in the inferred extrapolations of iF^rP to g^ = 0, where 
the G&S model gives results more consistent with the expectation ^,7(0) = 1. 

The presence of the p" in e^e~ — > vr^vr" is evident from the cross section measurements of 



DM2 1 28 1 near and above Mt-. For r lepton decay, the p" pole mass is near the endpoint of the 
M^-^o spectrum, thus making it difficult to observe. However, as with the p', its influence can be 
observed as an interference effect. While we obtain good fits without the p" meson, the x^ values 
for both K&S and G&S models are significantly improved when such a contribution is introduced. 
It can be seen in Fig. ^ that the p" fits also agree better with the data points in the 1.5-1.6 GeV 
region which were excluded from the fits. 
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FIG. 9. Alternate fits to the fully corrected M^-^o distribution in r^ -^ tt^tt^Vj. events (points). 
The solid curve overlaid represents the results of the fit to the K&S model with p" contribution 
included. The dashed curve represents the fit to the G&S model, also with the p" . The dotted 
curve is the fit to the K&S model, including the Blatt-Weisskopf barrier factor but no p" . The inset 
shows the low mass region where the differences between the models are most significant. Here we 
plot l-Fvrl^ (see Sec. VII B ), eliminating the purely kinematic factors which cause rapid variation in 
M^-^o in this region. 
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The relative phases of /?, p' and p" which we have obtained as (+ — +) are consistent with 
expectations from some models 1 1^,291. We also find that including the p" has a significant impact 



on our measurement of the p' parameters. In particular, values for the p' mass are closer to those 



based on other decay modes p7| , pO[ when the p" is included. 

We have also modified the energy dependence of the p and p' widths by including the Blatt- 
Weisskopf factor (see Eq. P) . The results for the nominal K&S fit function (with no p" contribution) 
modified in this way are given in Table || and shown as the dotted curve in Fig. p. Values obtained 
for the p and p' range parameters Rp and Rpi are consistent with expectations. This fit yields a 
smaller x^ per degree of freedom than other fits that also do not include the p(1700). However, 
this function as implemented does not yield a normalization consistent with -F7r(0) = 1. 

VI. SYSTEMATIC ERRORS 

Systematic errors are listed in Table |II[ These have been determined for the nominal K&S 
fit, however they are representative of those associated with G&S type fits as well. We have not 
included uncertainties associated with model dependence in our systematic error assessment. We 
discuss the most significant sources of error below. 

TABLE II. Systematic errors, in MeV except those for /3, which is dimensionless. 



Source 


Mp 


Tp 


/3 


Mp. 


Tp' 


Backgrounds 


0.1 


0.4 


0.002 


1 


5 


Bin Migration 


0.3 


0.5 


0.001 


5 


16 


Momentum Scale 


0.2 


0.1 


< 0.001 


1 


1 


Energy Scale 


0.8 


0.2 


0.002 


2 


6 


Acceptance 


0.2 


0.2 


0.004 


5 


15 


Fit Procedure 


0.1 


0.2 


< 0.001 


1 


1 


Total Syst. 


0.9 


0.7 


0.005 


8 


23 


Stat. Error 


0.5 


1.1 


0.007 


7 


26 



A. Background Subtraction 

We rescale the individual background components by amounts consistent with uncertainties 
in branching fractions and detection efficiency. Typically, this is ±(5-10)% of the nominal. The 
largest change in fit parameters comes about by varying the t~ -^ K~ti^Vt contribution in this 
way. 

B. Bin Migration Correction 

The resolution on M^-^o is dominated by the photon energy measurement. The excellent 
agreement between data and MC seen in the vr" mass spectra in Figure |l| gives us confidence that 
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this resolution function is adequately simulated. That the MC correctly models radiative effects 
is tested by comparing characteristics of photon-like showers accompanying the 7r~7r*^ system with 
those in the data. 

The method applied to unfold these effects relies on the approximate similarity between the 
M^-^o line shape used as input to the MC and the true one (see Appendix^. We have estimated 
the bias resulting from the inaccuracy of this approximation, and find it to be small. Considering 
this bias and possible errors in the modeling of the effects themselves, we conservatively arrive at 
the uncertainties shown in Table ||. For reference, failure to correct for bin migration results in 
values for Mp and Tp which are 2.6 MeV lower and 4.6 MeV higher, respectively. Thus we believe 
we understand this correction to ~ 10% of itself. 



C. Energy and Momentum Scales 

Although the vr'' mass resolution function is well modeled, even a small error in the absolute 
energy calibration of the calorimeter can have an effect on Mp without grossly disturbing the agree- 
ment in Figure 111 Based on studies of photons from various processes, we believe this calibration 



to be good to better than 0.3% |31]. Scaling the photon energies in the data by a factor of 1 ±0.003 
and fitting the M^-^o spectra thus obtained results in the uncertainties shown in Table O. 

The same procedure is used to assess the systematic error due to the absolute momentum scale 
uncertainty. Studies of D meson decays and e'^e~ -^ ^^~^ t^~ events have demonstrated that the 
momentum scale uncertainty is below 0.05% |32]. 



D. Other Sources of Error 

To estimate errors associated with the mass dependence of the detection efficiency, we modified 
the shape of the acceptance correction distribution (Fig. |^ by amounts suggested by the variation 
of selection criteria. We also performed the full analysis after varying cuts on 5-y-y, minimum photon 
energy, and the vr" — > 77 decay angle. Uncertainties associated with the fitting procedure were 
evaluated by performing fits to generator-level Monte Carlo spectra (no detector effects simulated). 
We also fit the uncorrected data and MC spectra, using the observed shifts in fit parameters for 
the MC sample to correct the parameters obtained from the data sample. This procedure yielded 
results that were in agreement with our nominal procedure. As a final cross-check, we split the 
data sample according to the tag decay. The results obtained for the three tags were in agreement 
with each other. 



VII. COMPARISON WITH DATA FROM OTHER EXPERIMENTS 

A. Comparison with Fits by Other Experiments 

The results from our fits can be compared with similar fits by ALEPH Q to their data, as 
well as with fits by other authors |13,|l4,l^ to various compilations of e^e~ — > 7r"*"7r~ data. For 



illustration purposes, some representative comparisons are given in Table III. For the e'^e~ data we 
give the results of the fits carried out by Kiihn and Santamaria |^] to data below ^/s ~ 1.6 GeV. 
These authors did not quote uncertainties on the parameters they obtained. However, based on 
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TABLE III. Comparison of results of fits using the K&S and G&S models [p" included with 
Mp//, Tp„ = 1700, 235 MeV) to (a) CLEO r data, (b) ALEPH r data § and (c) e+e" -^ vr+vr- 



data (denoted 'K&S Fit', from fits in Ref. JlSl). Units are MeV for all fit parameters, except (5, 
7 and |-F7r(0)P which are dimensionless. Shown in parentheses are the experimental uncertainties 
(statistical plus systematic). The parameter associated with p-co interference in the e^e~ data is 
not shown. 



Parameter 



CLEO 



K&S Model 
ALEPH 



K&S Fit 



CLEO 



G&S Model 
ALEPH 



K&S Fit 



Mp 



\F.m 



774.6(1.1) 


774.9 (0.9) 


773 


775.1(1.1) 


776.4(0.9) 


776 


149.0(1.4) 


144.2(1.5) 


144 


150.4(1.4) 


150.5(1.6) 


151 


-0.167(10) 


-0.094 (7) 


-0.103 


-0.121(10) 


-0.077(8) 


-0.052 


1408 (14) 


1363(15) 


1320 


1406(15) 


1400 (16) 


1330 


502 (39) 


= 310 


390 


455 (41) 


= 310 


270 


0.050(10) 


-0.015(8) 


-0.037 


0.032 (9) 


0.001 (9) 


-0.031 


1.14(2) 


= 1 


= 1 


1.03(2) 


= 1 


= 1 



xVdof 



23.2/23 



81/65 



136/132 



22.9/23 



54/65 



151/132 



similar fits performed by Barkov et al. [|18[ and by us, we believe these uncertainties to be similar 
in magnitude to those from the r decay data fits. For example, averaging over models, Ref. |18| 



obtains Mp = 775.9 ± 0.8 ± 0.8 GeV, Tp = 150.5 ± 1.6 ± 2.5 GeV, where the first error represents 
that due to statistics plus systematics, and the second error gives the uncertainty due to model 
dependence. 

The comparisons given in Table |T| are not meant to be rigorous. Due to the different assump- 
tions made by different authors, a systematic comparison is not possible. For example, the choice 
of whether to enforce the -F7r(0) = 1 condition for the K&S model has a strong impact on the fit 



parameters and the goodness of fit for the r data (see Sec. VB). In the G&S model this choice 
is less crucial since the values of |F7r(0)p obtained when it is allowed to float are closer to unity 
than they are in the K&S model. In addition, the fits to the e^e~ data shown do not include the 
high-^s data from DM2 [28|. With this data included, fits we have carried out yield results that 
are considerably different from those with this data excluded. The applicability of any given model 
across the full range of ^/s accessible to experiments has not been demonstrated. 

Finally, small differences between charged and neutral p meson parameters are expected, due 
to manifestations of isospin violation such as the 7r~/7r'^ mass difference. Refs. Q and ^J include 
discussions of this issue. Additional comments are given in Appendix [Q. 



With these caveats, Table 111 demonstrates general agreement between CLEO and ALEPH 
data, and between the r and e~^e~ data, supporting the applicability of CVC. In particular, the 
Gounaris-Sakurai fits show very good agreement between CLEO and ALEPH data. Within the 
K&S model, the smaller values of Tp obtained for the ALEPH and e"'"e~ data are likely a con- 
sequence of the i^7r(0) = 1 constraint, as similar values are also obtained when this constraint is 
imposed for the CLEO data (see Sec. [VB| ). 



B. Model Independent Comparisons 
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1. Comparison of CLEO and ALEPH data 

It is also useful to compare data from different experiments in a model-independent way. The 
ALEPH Collaboration has published their corrected M^_ o spectrum H, which may be compared 
directly with the corresponding CLEO spectrum. Given the differences in binning, a visually 
useful comparison is obtained by determining the bin-by-bin deviations of the two data sets from 
the prediction of a given model. For the prediction we use the results from the fit of the CLEO data 
to the G&S model, including p" . In Fig. KU, we plot the deviations of the CLEO and ALEPH data 
from this model as a function of M^-^o. The clustering of the CLEO points around zero is expected 
given the goodness of the fit. The ALEPH points also cluster around zero, although systematic 
deviations may be present at low and high masses. The significant correlations among the ALEPH 
points make this difficult to establish from the figure. Independent of possible deviations from the 
model, however, the CLEO and ALEPH points show agreement over the entire spectrum. 

2. Comparison of CLEO and e^e^ data 



Direct comparison, independent of model, of r and e^e~ data can also be made as a further 
test of CVC. To do this we represent the fully corrected CLEO vr^vr'^ mass spectrum in terms 
of I-Fjrp which can be compared with the e~^e~ data directly. Using the leptonic r decay width 
Fg = GpM^S']^]y/1927r^ and the ratio of r decay branching fractions to 'K~'K^Vr (B^^o) and e~Vf.VT 
{Bf.) for normalization, we derive the spectral function averaged over each mass bin from Eq. ||: 
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U-fM) = ^^ ^' ^^ ^ 1 -^ (15) 

^ '' B, 127r iKdP SY^ M, (M2 - MfY {M? + 2Mi) N AM,' ^ ' 

where Mj is the central value of M^-^o for the i*^ bin, and the quantity 1/N N^/AMi is the 
number of entries {N,) in the i^^ bin of the corrected measured 7r~7r'^ mass spectrum, divided by 
the total number of entries (A^) and the bin width (AMj). We use world average values p^ for 
the branching fractions, B^^o = (25.32 it 0.15) % and B^ = (17.81 it 0.07) %, for the r lepton mass 
Mr = 1777.05l|]|^ MeV, and for the CKM element |Kd| = 0.9752 ± 0.08. S1;y^, represents the 
radiative corrections to leptonic r decay, evaluated to be 0.996 (see Ref. 0). The overall radiative 
correction factor Sew = •S'^V/'^EVy ^^ estimated to be 1.0194 @, based mainly on the logarithmic 
terms associated with short-distance diagrams involving loops containing bosons, which differ for 
leptonic and hadronic final states. Non-logarithmic terms have not yet been computed for the tttt 
final state, but are expected to be small. 

\Ft^{M)\'^ is computed from v{M) using Eq. ^ A small (< 1%) correction is made to represent 
each point as a measurement at the central value of its mass bin, rather than as an average over 
the bin. This is done by employing the results from the fits described earlier to estimate the effect 
of the line shape variation across each of the histogram bins. Although this correction depends on 
the model used, this model-dependence is negligible relative to experimental errors. The factors 
appearing in Eqs. ^ and |l^ are also used to recast the normalization parameters determined in the 
fits described in Sec. ^ in terms of |-F,r(0)P, as shown in Table |. 

In Figure 0, the values of iF^rp derived from the corrected CLEO M^-^o spectrum are plotted 
along with the with e+e" -^ vr+vr- data from CMD ]l8|, CMD-2 |3|, DM1 H, DM2 |2§, 



OLYA 1 18], NA7 [^, as well as Adone |3£-38| experiments and other VEPP |18,3^ experiments. 
The r data follow the e'^e" data shape well, except in the region where p-oj interference affects the 
e^e~ data. However, the r data tend to lie above the e^e~ data throughout most of the range in 



^/s. This is illustrated in Fig. 12, where we plot the fractional difference between measured values 
of |-F7r('S)P from the e^e~ data and the prediction from the G&S fit (including p") to the CLEO 
r data. To obtain this prediction, we modified the G&S fit function for use with e~^e~ data. We 
allowed the normalization and a parameter associated with p-uj interference to float to determine 
the latter, and then readjusted the normalization to give that shown in Table |. 



The deviation from zero in Fig. 12 is consistent with the deviation from unity of the value 
of |F7r(0)p = 1.03 inferred from the fit to the CLEO r data. (Fitting the e+e~ data with the 
normalization floating yields values of |-F7r(0)P ~ 1.00.) It is also consistent with the observation 
made by Eidelman and Ivanchenko |40| that the measured r^ -^ tt^-k^Vt branching fraction is 
larger than that expected from application of CVC to the e^e~ — > vr+vr^ data by 3.2 it 1.4%. The 
consistency of these three indications of CVC violation is not accidental since they all involve (1) 
application of CVC to largely the same e'^e" data, as well as the use of (2) the world average 
value for the r decay branching fractions |^^ and (3) the same electroweak radiative correction 
factor [0,^] to normalize the r decay spectral function, as described above. Known sources of 
error in the first two of these components have been included in the ±1.4% error above. Recent 
estimates suggest that the uncertainty associated with additional radiative correction factors not 
yet computed could be as large as ±0.4% (see Refs. |^,p|). Deviations associated with isospin- 
violating effects are expected to be small, but carry an uncertainty not reflected above. 
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FIG. 11. Comparison of I-FttP as determined from CLEO-II r^ -^ tt^tt^Ut data (filled circles), 
with that obtained from e~^e~ -^ TT~^7r~ cross sections (other symbols) from CMD, CMD-2, DM1, 
DM2, OLYA, and NA7, as well as Adone and other VEPP experiments. The inset is a blow-up of 
the region near the p peak, where p-uj interference is evident in the e~^e~ data. 
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FIG. 12. Difference between l-F-n-p as determined from e^e~ — > vr^vr" data and that inferred 
from fit of the CLEO r~ — > tt'tt^u-j- data to the Gounaris-Sakurai model, with p" included, divided 
by the fit value. The different symbols represent the data points from the e~^e~ experiments 
operating in the p peak region. 
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VIII. IMPLICATIONS FOR THE MUON ANOMALOUS MAGNETIC MOMENT 

A. Introduction and Motivation 

The muon magnetic moment anomaly a^ = {g^ — 2)/2 has an experimental value of 
(11659 230 lb 84) x 10^^'^ [p7,42]. It receives significant contributions from uncalculable hadronic 



vacuum-polarization radiative corrections, estimated recently [^ to be a^'^ = (692.4ib6.2) x 10"^'^. 
The error on this contribution is the most significant source of uncertainty in the Standard Model 
prediction for a^. Experiment E821 at Brookhaven National Laboratory, currently under way, is 
aiming at a precision of ~ 4 x 10^^*^ in order to be sensitive to the weak interaction contribution 
at 15 X 10~^^, as well as possible contributions from new physics. Improved knowledge of a^'^ is 
important for interpreting experimental results on a^. In this section, we discuss the implications 
of our data on r~ -^ TT~7r^i'r for the determination of a?^ , following the first such treatment of r 
decay data, by Alemany, Davier and Hocker [Q who used ALEPH data Q. 

The value of a,, is related to the e'^e~ anihilation cross section to hadronic final states via 



the dispersion integral [44| 



had ^ ^errM ^ ^ ^(^^ ^(^^^ (^g^ 

'^ VT J4M2 S 

where v{s) denotes the inclusive hadronic spectral function, and aem,(0) is the fine structure con- 
stant. K{s) is the QED kernel 

K{s)=x^ I'l _ ^"J + (1 + ^)2 U + X^ {\n{l + x) - X -'^\ + {^-^^ x^lnx, (17) 

where x = (1 — /3^)/(l + /?^), with /3^ = (1 — 4M^/s)^/^. Consequently, predictions for a^"''^ are 
based primarily on e'^e" data. By virtue of the form of K{s) and the 1/s factor in the integrand 
above, measurements at low values of ^/s where the tt+tt" final state dominates are particularly 
significant. 



B. The Impact of r Decay Data 

The smallness of the present error on a^'^ reflects the use of CVC and r decay data from 
ALEPH Q to improve the precision relative to that obtained based on e^e~ data alone. Specifically, 
the authors of Ref. [^ determine the tt+tt" contribution over the interval ^/s = 0.320 — 2.125 GeV 
to be 



(495.86 ± 12.46) x 10"^° (e+e" data only) 
(500.81 ± 6.03) X 10-1° {e+e- + r data) 



ar[0.320, 2.125] = ^ in + — --^^ , (ig) 



where the e+e~ data did not include the CMD-2 data [33| which were not available at that time. 
The combining of e^e~ and r data was carried out by averaging data points from different exper- 
iments prior to integrating the data. An advantage of this procedure is that it properly weights 
the r data more heavily in regions where it is more precise than the e^e~ data and vice- versa. 
Accordingly, the improvement in precision indicated in Eq. ^ is greater than that obtained by 
determining a^^ separately for r and e+e^ data and averaging the results. 
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In assessing the impact of the CLEO 7r~TT^ spectral function on the value and precision of 
o^ , we perform the integration in Eq. ^ using our data alone. Although this procedure does not 
possess the benefit described above, future changes to the central values or errors of external factors 
such as B^^o or Sew can be propagated easily. However, we note that a careful determination of 
^had -^Q^ijj combine the data from different experiments in some fashion, for example along the 
lines of the approach described in Ref. 0]. 

Details of our determination of aJJ' , including the method and corrections applied to account 
for isospin violating effects are presented in Appendix ^. Here, we report our results for the same 



interval in ^/s as in Eq. 18. We obtain 



a 



[0.320, 2.125] = (513.1 ± 2.1 ± 3.0 ± 4.5) x 10"^", (19) 



where the first error is due to statistics, the second error is due to experimental systematic uncer- 
tainies, and the third error reflects the uncertainties in externally determined quantities entering 
the calculation. These sources of error are discussed in Appendix |C|. 



The deviation between the CLEO result above and the result from Ref. Q given in Eq. 18 
based on e^e~ — > vr+vr" data only is again indicative of disagreement in the normalizations, rather 
than in the shapes, of the e~^e~ and r spectral functions as discussed in the previous sections. 
While the magnitude of this deviation is not significant on the scale of the reported errors, both it 
and the errors are greater than the projected precision of the BNL E821 experiment. Continued 
efforts to precisely determine ayf" will be needed to help interpret the results of the Brookhaven 
experiment. 

IX. SUMMARY 

From a sample of approximately 87000 decays (after subtraction of backgrounds) of the type 
T" -^ TT^TT^Ur, wc havc investigated structure in the Tr^vr*^ invariant mass spectrum. From fits to 
several models we have obtained parameters and relative couplings for the p(770) and p(1450) res- 
onances. Within the Kiihn and Santamaria model ||I^ with no p(1700) contribution, the precisions 
on the p{770) mass and width are 1.0 and 1.4 MeV, respectively. These precisions are comparable 
to those obtained by ALEPH |^ using the same r decay mode, as well as those obtained from fits 
to low-energy e~^e~ -^ tt+tt" data p|, p^Jl^ . 

We find that a successful description of our M^-^o spectrum requires interfering Breit-Wigner 
line shapes associated with p(770) and p(1450) resonances. Fits including additional contributions 



from the /o(1700) resonance are preferred, consistent with the observations of DM2 |28| in e~^e~ -^ 
vr+vr^. Fits of comparable quality are also obtained without a p" contribution when the mass- 
dependence of the p and p' widths is modified to include Blatt-Weisskopf barrier factors. We find 
that Gounaris-Sakurai [|ll| type fits yield extrapolations to g^ = that are higher than but roughly 
consistent with -F^(O) = 1. The Kiihn-Santamaria [^] type fits lead to significantly higher values, 
such that imposing ^,^(0) = 1 can lead to biased values for p and p' parameters. 

Quantitatively, the central values for the precisely determined p(770) parameters are difficult 
to compare with those from fits by others. This is due to severe model dependence and the strong 
influence of data far from the p peak. Qualitatively, the shape of our M^-^o spectral function 
agrees well with that obtained by ALEPH ^. It also agrees well with the e"'"e~ -^ tt+tt" data, 
supporting the applicability of CVC. Using our spectral function, we have employed CVC to infer 
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the significant component of the hadronic contribution to the muon g — 2 factor associated with 
the e+e" -^ vr+vr" cross section, obtaining a^''[0.320, 2.125] = (513.1 ± 5.8) x 10"^°. 

However, we also observe indications of discrepancies between the overall normalization of r 



and e^e~ data, as pointed out by Eidelman and Ivanchenko |4C]. These appear in the model 
independent comparison of our spectrum with the e~^e~ cross section measurements, as well as 
in the values for |F7r(0)p inferred from fits to our spectrum which are larger than unity, and in 
our high value for aj^^. Though larger than the deviations expected due to known sources of 
isospin violation, these could also arise from experimental errors in the r decay branching fraction 
measurements or the normalization of the e~^e~ — > vr^vr" cross section measurements, or from 
theoretical uncertainties in the estimates of radiative corrections, or from all of these sources. At 
present the deviations are not significant on the scale of the reported errors. We hope that new 
data from Novosibirsk, BEPC, and the B-factory (PEP-II, KEK-B and CESR-III) storage rings 
will shed light on this issue in the near future. 
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APPENDIX A: UNFOLDING OF RESOLUTION AND RADIATIVE EFFECTS 

The unfolding procedure is used to correct for bin migration in the observed M^^^o spectrum 
due to resolution smearing and radiative effects. In this Appendix, we describe the procedure used 
in this analysis. In the discussion that follows, matrices are denoted in boldface, while non-bold 
symbols represent vectors with the contents of binned histograms. 

The unfolding procedure makes use of the Monte Carlo to characterize the bin migration. One 
may construct a migration matrix P which gives the probability that an event generated with a 
given 7r^7r'^(7) mass is reconstructed with a given tt^tt^ mass: 

^MC^p^MC^ (Al) 

where G^^ represents the generated 7r~7r°(7) mass spectrum and B}^'~" represents the reconstructed 
one. It is possible to apply the inverse of P to the spectrum observed in the data i?"^*^ to obtain 
an unfolded spectrum [7°^*^ that provides an estimate for the parent spectrum. However, such a 
procedure is not robust with respect to statistical fluctuations entering the determination of P, 
and can yield spectra with unphysically large point-by-point fluctuations. 
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The corrected ALEPH spectrum |^] was derived using a method based on singular value de- 
composition of the migration matrix [^ to mitigate this effect .Q Here, we use an iterative method 
that relies on the smallness of the bin migrations and the approximate similarity between the 
reconstructed spectra from the data {R'^'^^'^) and Monte Carlo {RJ^'^) samples. We construct the 
matrix P', which gives the fraction of MC events with a given reconstructed tt~it^ mass that were 
generated with a given 7r~7r''(7) mass. With this matrix, G = P' R is satisfied, but P' is 
not equal to P~^. Successive application of P' to the observed data spectrum gives an estimate 
for the parent distribution, according to: 



oo 



rydata _ r>/ odata _j_ V^ (i _ pi p\k pi / ndata _ E>MC\ ( A2'\ 



For small bin migration probabilities, all elements of the matrix (1 — P'P) are small: these 
elements are the 'expansion parameters' in the series above. A simpler form for Eq. ^^ is obtained 
by recognizing that the quantity {1 — P' P)P' R is identically zero. As written, however, Eq. A2 



illustrates that when the series is truncated, deviation from the results of the full expansion vanishes 
as (i?^^*^ - R^^) -^ 0. 

With the similarity between the observed data and MC spectra plotted in Fig. ||, we find it 
sufficient to ignore terms with A; > 1 in Eq. E^. The k = 1 term has a noticeable effect on several 



of the fit parameters reported in Section VB , In particular, ignoring this term leads to a value for 
Tp that is 0.5 MeV smaller than that from our nominal fit. Higher order terms have no significant 
impact on any of the parameters. 

APPENDIX B: THE CORRECTED 77-77° MASS SPECTRUM 
In this appendix we present the fully corrected CLEO M^^o spectrum in tabular form. The 



spectrum is given as a compilation of event yields for each mass bin in Table IV, normalized so that 
the sum of entries over all bins is unity. Also given are the square roots of the diagonal elements 
of the covariance matrix (statistical errors only) . 



The statistical errors given in Table IV do not reflect the correlations between data points that 
were introduced by the bin migration correction. Although the correlations are not large, proper 
treatment of these data necessitates use of the covariance matrix. The full covariance matrix y is a 
43 X 43 symmetric matrix. In Table 0, we present the correlation coefficients pij = Vij/iViiVjjY''^ 
for bins i and j, i > j, for i, j where pij > 0.0015. The coefficients shown are the statistical 



correlations only. Both the corrected spectrum as given in Table [V and the full covariance matrix 
are available electronically |23]. 



APPENDIX C: DETERMINATION OF A^^" FROM CLEO DATA 

In this appendix, we present the details of our determination of aJJ^ over the range ^/s = 0.320 
to 2.125 GeV. Quantities entering this determination are summarized in Table VI. We describe the 



"^For fits to models, ALEPH performed fits of their uncorrected spectrum to convolutions of the 
experimental effects with the functions describing the models, avoiding this problem altogether. 
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TABLE IV. The corrected CLEO M^^o spectrum given in terms of event yield as a function of 
mass bin. The spectrum is normahzed so that the number of entries sums to unity. The numbers in 
parentheses denote the square roots of the diagonal elements of the covariance matrix (statistical 
errors only). Note that the number of entries jumps at the 1.00-1.05 GeV bin, where the bin size 
is increased from 0.025 to 0.050 GeV. 



Bin 


Mass Range 


Entries 


Bin 


Mass Range 


Entries 


No. 


(GeV) 


(10-4) 


No. 


(GeV) 


(10-4) 


1 


0.275-0.300 


1.5 (1.4) 








2 


0.300-0.325 


8.0 (2.5) 








3 


0.325-0.350 


6.0 (2.6) 








4 


0.350-0.375 


8.5 (2.3) 


24 


0.850-0.875 


446.7 (6.6) 


5 


0.375-0.400 


15.6 (2.6) 


25 


0.875-0.900 


326.2 (5.5) 


6 


0.400-0.425 


16.2 (2.9) 


26 


0.900-0.925 


262.1 (4.9) 


7 


0.425-0.450 


24.9 (3.1) 


27 


0.925-0.950 


207.3 (4.3) 


8 


0.450-0.475 


41.4 (3.4) 


28 


0.950-0.975 


158.8 (3.7) 


9 


0.475-0.500 


50.6 (3.7) 


29 


0.975-1.000 


129.6 (3.4) 


10 


0.500-0.525 


60.9 (4.0) 


30 


1.000-1.050 


202.8 (4.8) 


11 


0.525-0.550 


79.8 (4.4) 


31 


1.050-1.100 


151.0 (4.1) 


12 


0.550-0.575 


107.4 (4.7) 


32 


1.100-1.150 


111.0 (3.4) 


13 


0.575-0.600 


144.3 (5.2) 


33 


1.150-1.200 


87.0 (3.0) 


14 


0.600-0.625 


204.5 (5.9) 


34 


1.200-1.250 


63.9 (2.5) 


15 


0.625-0.650 


269.1 (6.5) 


35 


1.250-1.300 


42.7 (2.0) 


16 


0.650-0.675 


385.8 (7.5) 


36 


1.300-1.350 


29.2 (1.7) 


17 


0.675-0.700 


571.5 (8.7) 


37 


1.350-1.400 


18.1 (1.3) 


18 


0.700-0.725 


826.8 (10.1) 


38 


1.400-1.450 


6.98 (0.84) 


19 


0.725-0.750 


1078.4 (11.3) 


39 


1.450-1.500 


2.91 (0.59) 


20 


0.750-0.775 


1228.1 (11.8) 


40 


1.500-1.550 


0.71 (0.32) 


21 


0.775-0.800 


1114.7 (11.0) 


41 


1.550-1.600 


0.59 (0.25) 


22 


0.800-0.825 


878.1 (9.6) 


42 


1.600-1.650 


0.68 (0.26) 


23 


0.825-0.850 


629.3 (7.9) 


43 


1.650-1.700 


0.28 (0.21) 
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TABLE V. 


Correlation coefficients pij for contents of bins i 


and J 


(see Table |vl) of the fully 


corrected CLEO M 


ttttO 


spectrum 


, for 


i>: 


'. Coefficients with values less than 0.0015 


are 


not shown. 


These coefficients refiect statistical correlations 


only. 


















i 3 


P^j 


i 


j 


Pij 


i 


3 


Pij 


i 


3 


Pi3 


i 


3 


Pi3 


i 


3 


Pij 


2 1 


-.002 


10 


8 


-.057 


16 


14 - 


-.084 


21 


19 


-.116 


28 


27 


.299 


38 


35 


-.011 


3 1 


-.031 


10 


9 


.085 


16 


15 


.183 


21 


20 


.205 


29 


25 


-.006 


38 


36 


-.074 


3 2 


.042 


11 


6 


-.002 


17 


3 


.002 


22 


18 


-.003 


29 


26 


-.046 


38 


37 


.114 


4 1 


.003 


11 


8 


-.013 


17 


11 - 


-.002 


22 


19 


-.033 


29 


27 


-.102 


39 


35 


.002 


4 2 


-.025 


11 


9 


-.064 


17 


12 - 


-.003 


22 


20 


-.113 


29 


28 


.317 


39 


36 


-.009 


4 3 


.053 


11 


10 


.120 


17 


13 - 


-.003 


22 


21 


.219 


30 


26 


-.004 


39 


37 


-.082 


5 1 


.007 


12 


2 


.002 


17 


14 - 


-.023 


23 


1 


.003 


30 


27 


-.034 


39 


38 


.103 


5 2 


-.008 


12 


7 


-.003 


17 


15 - 


-.090 


23 


19 


-.003 


30 


28 


-.094 


40 


37 


-.015 


5 3 


-.034 


12 


9 


-.014 


17 


16 


.201 


23 


20 


-.036 


30 


29 


.123 


40 


38 


-.079 


5 4 


.012 


12 


10 


-.076 


18 


13 - 


-.002 


23 


21 


-.107 


31 


28 


-.019 


40 


39 


.166 


6 1 


-.004 


12 


11 


.145 


18 


14 - 


-.005 


23 


22 


.241 


31 


29 


-.071 


41 


38 


-.011 


6 2 


.002 


13 


9 


-.002 


18 


15 - 


-.025 


24 


20 


-.003 


31 


30 


.029 


41 


39 


-.090 


6 3 


-.005 


13 


10 


-.017 


18 


16 - 


-.096 


24 


21 


-.037 


32 


29 


-.007 


41 


40 


.152 


6 4 


-.031 


13 


11 


-.075 


18 


17 


.206 


24 


22 


-.104 


32 


30 


-.061 


42 


36 


-.002 


6 5 


.027 


13 


12 


.150 


19 


12 - 


-.002 


24 


23 


.255 


32 


31 


.046 


42 


39 


-.013 


7 4 


-.006 


14 


3 


.002 


19 


14 - 


-.002 


25 


21 


-.003 


33 


30 


-.005 


42 


40 


-.094 


7 5 


-.041 


14 


9 


-.002 


19 


15 - 


-.004 


25 


22 


-.039 


33 


31 


-.066 


42 


41 


.226 


7 6 


.059 


14 


10 


-.002 


19 


16 - 


-.028 


25 


23 


-.102 


33 


32 


.056 


43 


18 


.002 


8 2 


-.002 


14 


11 


-.019 


19 


17 - 


-.103 


25 


24 


.261 


34 


31 


-.006 


43 


30 


-.002 


8 4 


-.002 


14 


12 


-.078 


19 


18 


.207 


26 


22 


-.004 


34 


32 


-.071 


43 


31 


.006 


8 5 


-.009 


14 


13 


.155 


20 


13 - 


-.002 


26 


23 


-.041 


34 


33 


.068 


43 


35 


-.003 


8 6 


-.053 


15 


11 


-.003 


20 


14 - 


-.002 


26 


24 


-.104 


35 


32 


-.007 


43 


36 


.007 


8 7 


.067 


15 


12 


-.019 


20 


15 - 


-.002 


26 


25 


.270 


35 


33 


-.073 


43 


37 


-.002 


9 3 


-.003 


15 


13 


-.082 


20 


16 - 


-.003 


27 


23 


-.004 


35 


34 


.074 


43 


40 


-.023 


9 4 


-.003 


15 


14 


.178 


20 


17 - 


-.030 


27 


24 


-.042 


36 


33 


-.009 


43 


41 


-.111 


9 6 


-.010 


16 


2 


.003 


20 


18 - 


-.111 


27 


25 


-.106 


36 


34 


-.076 


43 


42 


.246 


9 7 


-.056 


16 


3 


-.002 


20 


19 


.200 


27 


26 


.276 


36 


35 


.076 








9 8 


.075 


16 


10 


-.002 


21 


16 - 


-.002 


28 


24 


-.004 


37 


34 


-.008 








10 5 


-.002 


16 


12 


-.003 


21 


17 - 


-.003 


28 


25 


-.044 


37 


35 


-.073 








10 7 


-.011 


16 


13 


-.022 


21 


18 - 


-.032 


28 


26 


-.103 


37 


36 


.081 
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integration procedure, corrections for isospin-violating effects, and the evaluation of errors below. 
Finally we discuss the implications of our results. 



TABLE VI. Components of the derivation of a^'^ [0.320, 2.125] from the CLEO r" 
hadronic mass spectrum. 



vr tt^Vt 



Quantity 



Value 



Correction to 



Integration of v{s) (Eg. |16D 
g;^^ [0.320, 2.125] (raw) 



Uncertainty on 



(514.8 ±2.1) X 10" 



-10 



±2.1 



Normalization Factors 

B{t- -^ TT-TT^Vr) 

Sew [P] 



VudM 



(25.32 ±0.15)% 
(17.81 ± 0.07) % 
1.0194 ± 0.0040 
0.9752 ± 0.0008 



±3.0 
±2.0 
±2.0 
±0.8 



Correction Factors 
a {p-oj interference) 
AF. 



M^ 



M^ 



(1.71 ±0.06 ±0.20) X 10-3 

(±0.28 ± 0.39) % 

4.6 MeV 



±3.6 
±1.1 
-6.5 



±0.4 
±1.6 



Other Sources of Systematic Error 

Backgrounds 

Bin Migration 

Energy Scale 

Acceptance 

Integration Procedure 



±0.c 



±2.4 
±1.4 
±1.0 
±0.5 
±1.0 



1. Integration Procedure 



To evaluate aj^^ as given by Eq. ^, we perform a numerical integration employing the Gounaris- 
Sakurai model, with p" included, using the best fit parameters given in Table |. The externally 
measured quantities used to infer v'^^ from our corrected tt^tt^ mass spectrum are listed in Ta- 



ble VI. From this procedure we obtain, prior to application of the corrections described below. 



'[0.320, 2.125] = (514.8 ±2.1) x 10 



-10 



(CI) 



where the error is statistical only. 

The advantages of this procedure relative to direct integration of the data points (summing 
the histogram entries) are (1) mitigation of the effects of statistical fluctuations particularly in the 
low-mass bins and (2) operational simplicity. The disadvantages include possible biases associated 
with choice of model. We have checked this method of integrating the functional form of v^'^{s) 
by reproducing the CMD-2 evaluation [33] of a!^'^[0.61, 0.96] to within 0.6 x lO"^*^ of the value 
they obtained via direct integration of their data. We have also verified with our data that direct 
integration gives comparable results to those obtained by integrating the fit function. 
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2. Corrections 

Several corrections are needed to account for sources of isospin violation that bias the naive 
application of CVC. The corrections can be classified according to three quantities: (1) the magni- 
tude of the fhu interference arising from the isospin- violating electromagnetic decay uj — > vr+vr^; (2) 
possible isospin splittings between charged and neutral p meson masses and widths; and (3) kine- 



matic effects associated with the tt /tt^ mass difference. These corrections are listed in Table Vl, 
and are described below. 



a. Contributions from uj -^ vr^vr 
To account for the absence of p-oj interference in r data, we modify the G&S function to include 



it, introducing the parameter a (following the notation of Ref. |13|) to quantify the u; admixture, 
analogous to the parameters (3 and 7 which quantify the p' and p" amplitudes relative to that of 
the p meson. From fits to the e^e~ data, we find a = (1.71 it 0.06 it 0.20) x 10~^. Modifying our 
fit function in this way leads to an increase in a^^ by 3.6 x lO"^'^ units. 

b. p meson isospin splittings 

Following the authors of Ref. |Q , the charged-neutral p mass splitting, expected to be small, is 
taken to be zero. We use their evaluation |Q] of the pole width splitting AF^ = iT p- — rpo)/Fp = 
(2.8 lb 3.9) X 10"'^. The dominant source of this splitting is the 7r~/7r'^ mass difference, which gives 
rise to different kinematic factors for charged and neutral p decay (see Eq. g). Additional small 
differences between charged and neutral p meson decay also affect the widths, some of which are 
included in the above estimate. One possibly significant difference, not accounted for here, is the 
effect of final state Coulomb interactions in p^ -^ -k'^tt^' decay which would tend to enhance the p^ 
width 1^^. Finally, as in Ref. Q, we assume that the charged and neutral p' and p" parameters 
are the same. 

In determining the effect of the estimated AFp on a"^ , we modify Fp as it appears in the 
denominator of the expression for v'^'^{q^) (see, for example, Eq. 0). In the context of the K&S 
and G&S models, where the -F7r(0) = 1 constraint is enforced, the Fp factor does not appear 
explicitly in the numerator of the squared Breit-Wigner formula, unlike the general form for v'^'^{q^) 
given by Eq. |^. Modifying Fp in the denominator only leads to an additive correction to a^^ of 
(+1.1 lb 1.6) X 10^^^. This is contrary to the result of Ref. [^, in which the correction is given as 
-1.4 X 10-1°. 



c. Additional kinematic factors 

As mentioned above, the vr^/vr'^ mass difference contributes to AFp through the P-wave phase 
space factor (2^7^/^/^^)^ appearing in Eq. ^. This factor also characterizes the q^ dependence of the 
p width and affects the numerator as well as the denominator of v'^'"{q^) in the models considered. 
This effect influences the spectral function strongly at low values of q^, since that is where the 
values of p^r for charged and neutral p meson decay differ most significantly. Accounting for this 
difference leads to a decrease in a"^ by 6.5 x 10 
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3. Errors 



In this section, we discuss the sources of error indicated in Table Vl . 



a. Statistical Errors 

Since we use external measurements to normalize our spectral function, the statistical errors 
considered here are those associated with the bin-by-bin fluctuations in our mass spectrum. Statis- 
tical errors associated with the Monte Carlo based corrections for backgrounds, bin migration and 
acceptance also enter. We assess the overall statistical error by generating a large number of G&S 
parameter sets, with the parameters determined randomly about the central values returned by our 
nominal fit, weighted according to the covariance matrix returned by the fit, assuming Gaussian 
errors. We determine aV^ separately for each parameter set. The r.m.s. of the distribution of 
values was found to be 2.1 x 10"^'^. 



b. Internal Systematic Errors 

Internal systematic errors are those associated with our analysis of t~ -^ n^TT^i'r decays. They 
originate from the sources indicated in Sec. VI in the context of our fits to models of the TT~Tr^ 



mass spectrum. As with the statistical error, these errors pertain to the shape, rather than the 
normalization, of the spectral function. 

As expected from Table ||, the dominant sources are uncertainties associated with the back- 
ground subtraction and bin migration corrections. Possible biases in these corrections would tend 
to affect the low end and p peak regions of the mass spectrum, on which a'^'^ depends most sensi- 
tively. We have also considered energy scale and acceptance uncertainties. We have estimated the 



uncertainties associated with these sources, shown in Table VI, in the same ways as described in 



Sec. VI, 



We have also estimated the bias associated with the model dependence of the approach used to 
compute a^'^. This has been done by comparing values of a'^'^ obtained with different models, as 
well as by directly integrating the data points. We estimate an uncertainty of ±1.0 x 10^^*^ from 
this source. Adding this in quadrature with the errors described above yields an overall internal 
systematic error of ±3.0 x 10"^''. 



c. External Systematic Errors 

External systematic errors are those associated with the parameters used to infer aJJ^ from 
our corrected tttt^ mass spectrum. They include uncertainties associated with normalization fac- 
tors, of which B^^o, Be, Sew, and Vud contribute the dominant errors. They also include the 
uncertainties associated with the corrections for isospin-violating effects described in the previ- 



ous section. Adding the errors listed in Table VI for these sources in quadrature gives an overall 
external systematic error of ±4.5 x 10"^''. 
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4. Results and Discussion 



With the normahzation and correction factors listed in Table Vl, we obtain 
a;^''[0.320, 2.125] = (513.1 ± 2.1 ± 3.0 ± 4.5) x 10"^° 

where the first error is the statistical error, the second is the internal systematic error, and the third 
is the external systematic error. In the above expression, we have made explicit the dependence 
on the external normalization factors so as to facilitate incorporation of future measurements of 
these quantities. 

Since this evaluation of aj^ is independent of the e^e~ only estimate from Ref. [Q], we can 
perform a weighted average of the two results. From this, we obtain 

a;^''[0.320, 2.125] = (510.0 ± 5.3) x 10"^°. (C3) 

This determination does not include the ALEPH data [0], nor does it include the CMD-2 data j33| . 
This is larger than the value of (500.81 it 6.03) x lO"^'^ obtained by the authors of Ref. |Q], despite 
the apparent agreement of the CLEO and ALEPH 7r~7r'' mass spectra. This reflects in part the 
difference in the procedures of combining the data, as described above. 
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